Giant Anisotropy of Spin-Orbit Splitting at the Bismuth Surface 
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We investigate the bismuth (111) surface by means of time and angle resolved photoelectron 
spectroscopy. The parallel detection of the surface states below and above the Fermi level reveals 
a giant anisotropy of the Spin-Orbit (SO) spitting. These strong deviations from the Rashba-like 
coupling cannot be treated in k • p perturbation theory. Instead, first principle calculations could 
accurately reproduce the experimental dispersion of the electronic states. Our analysis shows that 
the giant anisotropy of the SO splitting is due to a large out-of plane buckling of the spin and orbital 
texture. 
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The realization of transistors via the transport of spin 
polarized electrons attracts the interest of the solid state 
community since 20 years [TJ[2]. In such devices, an ap- 
plied gate voltage induces a spin torque of the injected 
electrons via Spin-Orbit (SO) interaction. The energy 
scale of this effect is typically 1-10 meV in the semicon- 
ductor heterostructures [3], but reaches values 10 times 
larger at surfaces of systems containing heavy elements 
[IHSIIH]. Therefore, the latter are considered as valuable 
models for future spintronic applications [9]. This cross 
fertilizing field has been recently enriched by the discov- 
ery of protected edge states in topological insulators [10] . 

The effects of SO coupling at the surface of solid states 
materials have been reviewed by several authors [TTHT3] . 
As originally noticed by Rashba, the spin degeneracy of 
the electronic states is lifted by the breakdown of in- 
version symmetry [14]. Being a relativistic effect, the 
SO splitting arises from the asymmetry of the electronic 
wavefunction in proximity of the ionic cores [15] . Despite 
the complexity of this problem, the k • p perturbation 
theory provides the leading terms of the SO Hamiltonian 
for small electronic wavevectors. When the surface has 
Cs v symmetry, the first order term is indeed the Rashba 
hamiltonian aR(k y ). This interaction term 

leads to an isotropic spin splitting and chiral spin tex- 
ture. It reproduces correctly the Shockley states at the 
(111) surface of gold as well as the electronic properties of 
several heterostructures [4] [14] [16]. Nonetheless, higher 
order expansions become necessary in systems where the 
surface state cannot be modeled within the framework 
of the nearly-free-electron approximation. As an exam- 
ple, the topological insulators Bi2Te3 and Bi2Se3 show a 
trigonal warping which is linked to the third order term 
ftx ~ ky) a z [T7HT9] . As a consequence, the spin polariza- 
tion acquires an out-of-plane component [20] that would 
affect the spin transport in non-ballistic devices [21]. 

In this letter we show that the large deviations from 
the Rashba hamiltonian induce a giant anisotropy of the 
SO splitting at the (111) Bismuth surface. We make 
use of an ultrashort laser pulse to efficiently populate 



electronic states up to 0.5 eV above the Fermi Level. 
By these means, we observe that the SO splitting is 
ArM = 150 ± 10 meV along the T-M direction and in- 
creases by 300% for a rotation of 6° around the surface 
normal. We checked the SO anisotropy by ab initio cal- 
culations of the band structure, finding excellent agree- 
ment with the observed dispersion of electronic states. 
The analysis of the resulting wavefunctions indicates that 
the SO anisotropy is due to the buckling of the spin ori- 
entation out of the surface plane. By projecting the 
spin-polarized wave functions on a local basis, we find 
an out-of-plane component larger than 30%. Our result 
differs from previous DFT calculations [22 whereas cor- 
roborates the spin resolved measurements of the occupied 
electronic states [23]. In order to test the entanglement 
between the spin and orbital degrees of freedom [24], we 
also calculated the projections of the Kohn-Sham wave- 
functions on the atomic orbital set. It follows that spin 
and orbitals have an opposite and nearly proportional 
polarization. 

Photoelectron spectra have been collected with a pho- 
ton energy of 18 eV, energy resolution of 10 meV and 
angular resolution of 0.2 degrees. The bismuth (111) 
surface has been prepared by sputtering-annealing cycles 
of a polished monocrystal. Time resolved photoelectron 
spectroscopy experiments have been performed with the 
FemtoARPES setup, using a Ti: Sapphire laser that gen- 
erates 35 fs pulses centered at 790 nm with repetition 
rate of 250 kHz. Part of the beam is employed to gener- 
ate the fourth harmonic by a cascade of frequency mixing 
in BBO crystals (/3-BaB 2 4 ) [25 . The 197.5 nm probe 
and the 790 nm pump are focused on the sample with a 
spot diameter of 100 /im and 200 /im, respectively. Their 
cross-correlation in a BBO crystal has a full width at half 
maximum (FWHM) of 80 fs. The bandwidth of the 197.5 
nm beam (6.3 eV) limits the overall energy resolution of 
TRPES spectra to 60 meV. All the time resolved mea- 
surements have been performed with incident pumping 
fluence of 0.6 mJ/cm 2 . 

The DFT calculations have been done using the 
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FIG. 1: Photoelectron intensity map integrated in an energy 
window of 20 meV around the Fermi level. The data have 
been acquired with circularly polarized photons centered at 
18 eV. The red and blue line stand for the wavevectors cut of 
Fig. [3p,D and Fig. [4]\,B, respectively. 



"APW + lo " method of the WIEN2k code [26 with the 
SOI taken into account and the ex-corr. functional based 
on GGA [27]. The surface was modelled by a symmetric 
slab of 20 Bi layers repeated along [111] with a periodic 
gap of 10 A [5j[6]. The interlayer distances in the slab 
were set to the values determined by the X-ray diffraction 
tecnique [28 . 

Figure [I] shows a photoelectron intensity map acquired 
with 18 eV photons in an energy window of 20 meV cen- 
tered at the Fermi level. The crossing points of surface 
states generate an internal pocket of hexagonal shape and 
6 elongated lobes along the T-M direction. The resulting 
map is in excellent agreement with the original measure- 
ments of Ast et al. [29 and bears no resemblance to 
the Fermi surface of nearly- free-electron systems [4]. In 
a following work, Korotheev et al have reproduced the 
experimental band structure by performing first principle 
calculations that accounted for the SO coupling [5 , 6]. By 
these means, they proved that SO splitting has a major 
effect on the electronic properties [7]. 

In agreement with these results Fig. [2] shows that the 
SO interaction removes the Kramers degeneracy of the 
surface states other than the T point. We recall that 
these wavefunctions hybridize with bulk states in the 
near proximity of the zone center while they are truly 
evanescent for larger wavevectors [6] . The data show that 
both spin-polarized bands cross the Fermi level along the 
T-M direction (panel A) whereas a single band crosses 
Ep along T-K (panel B). Independently on the crystal- 
lographic axis, the higher lying band is above the Fermi 
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FIG. 2: A-B): Photoelectron intensity map acquired along 
the r-M direction (Panel A) and T-K direction (Panel B) is 
compared to the calculated band structure. The data have 
been acquired with circularly polarized photons centered at 
18 eV. The solid lines stand for the two surface states whereas 
the dashed line marks the upper boundary of the projected 
bulk states that are below Ef- 



level for parallel wavevector 0.07 < k\\ < 0.5 A -1 . Since 
these states are not accessible by standard ARPES ex- 
periments, the value of the SO splitting could be inferred 
only by the DFT calculations. 

In order to overcome this limitation, we measure the 
SO splitting by means of time resolved Photoelectron 
Spectroscopy. The bismuth surface is excited by an ultra- 
fast pulse centered at 1.6 eV and is subsequently probed 
by a delayed 6.3 eV pulse. Figure [3]A. shows the photo- 
electron intensity map acquired as a function of pump- 
probe delay for k\\ = 0.35 A -1 along the F-M direc- 
tion. The transient population of electronic states de- 
cays via emission of phonons on a timescale of 6 picosec- 
onds. Since the electron gas thermalizes within 100 fs, 
the occupation factor is approximately described by a 
Fermi-Dirac distribution. Figure [3^ displays the tempo- 
ral evolution of the electronic temperature after photoex- 
citation. More details on the procedure to extract this 
parameter can be found in the supplementary material 
of Ref. [30 . Notice in Fig. ^ that the SO-splitting 
does not depends appreciably on the pump-probe delay. 
We conclude that the potential barrier at the surface 
of bismuth is barely affected by the elevated electronic 
temperature. This finding is somehow unfortunate, as it 
precludes the possibility to modulate the SO coupling by 
light fields. Nonetheless, a spin dynamics may still be vis- 
ible on the doped interface of Bi2Se3 [31] . In this case, the 
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FIG. 3: A): Photoelectron intensity map as a function of 
pump-probe delay for k\\ — 0.35 A -1 along the T-M direc- 
tion. B): Temporal evolution of the effective electronic tem- 
perature. C-D): Photoelectron intensity map acquired along 
the T-M direction (red line in Fig. [1) just before (C) and 200 
fs after (D) the arrival of the pump beam. Solid lines stand 
for the calculated dispersion of the surface states while no 
projected bulk bands exist in between the two dashed lines. 
The data of this figure have been generated by ultrafast pulses 
of linearly polarized photons centred at 6.3 eV. The electric 
field polarization was nearly parallel to the surface plane and 
formed an angle a — 20° with respect to the Gamma-M di- 
rection. 



200- 



> 

CD 

E 

CD 



-200- 



-400- 




4 -4 
Azimuthal angle (Degrees) 







D 








Band 2 




0.6 












0.4 








0.2 








0.0 








0.2 




"P o 


-3.5L Z 


0.4 


T 1 


1 1 1 


I i 


-4 





4 



Azimuthal angle (Degrees) 



FIG. 4: A-B): Photoelectron intensity map as a function of 
azimuthal angle for k\\ — 0.35 A -1 (blue line in Fig. [I]) and 
pump-probe delay of 200 fs. The probing photons are po- 
larized either parallel (Panel A) or orthogonal (Panel B) to 
the T-M direction. The white lines stand for the calculated 
dispersion of the surface states. C-D): In-plane component of 
the spin (a y , red dashed line), out-of-plane component of the 
spin (<7 Z , solid blue line) and of the orbital momentum (L z , 
open circles) for the lower band (Panel C) and upper band 
(Panel D) as a function of azimuthal angle. The orbital mo- 
mentum L z has been multiplied by -3.5 for better comparison 
with the out-of-plane spin. 



electronic structure of quantum well states could change 
upon the photoinduced reduction of band bending. 

Figure [3p and [3p shows the photoelectron intensity 
maps aquired along the T-M direction for a pump-probe 
delay r = —200 fs and +200 fs, respectively. The ob- 
served SO splitting is almost constant in the measured 
fc|| interval. Our experimental estimate Ap^ = 140 ± 10 
meV is consistent with A^ M = 130 meV extracted from 
DFT calculations. The SO spitting is strongly reduced 
with respect to the atomic value but it is large for a de- 
vice that would operate at room temperature. Notice in 
Figure [3]A and [3p that the pump beam also populates 
a bulk derived band above the surface state. This struc- 
ture at 370 meV nearly matches the lower boundary of 
the projected bulk bands that are above Ep. 

In the following we focus on the giant anisotropy of the 
spin-orbit splitting. We show in Fig. [4] the photoelectron 
intensity maps acquired at feii = 0.35 A -1 and r = 200 



fs as a function of azimuthal angle ip. The probe pulse 
polarization was either parallel (panel A) or perpendicu- 
lar (panel B) to the T-M direction. Notice that the lower 
and higher lying bands are more intense for parallel and 
perpendicular polarization, respectively. We evince that 
the surface states have opposite parity with respect to the 
T-M plane. According to the selection rules of the dipole 
matrix elements the band below Ep is even while the 
upper one is odd. The splitting between the two bands 
is minimal along the T-M direction but strongly depends 
on the azimuthal angle. When cp reaches ±6° the SO 
splitting attains 500 meV, thus increasing by 300% with 
respect to Ap^. The DFT calculations display a simi- 
lar trend but slightly underestimate the splitting value. 
Such large anisotropy is due to: a) the very large SO 
coupling of the atomic bismuth and b) the existence of 
surface states with large electronic wavevectors [20]. As 
a term of comparison, the second condition is not verified 
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in Bi 2 Te3 and Bi 2 Se3, where the SO anisotropy induces 
only a minor warping of the Dirac cone [T7HT9] . 

The question arises as which spin texture originates 
from to the measured band structure. On the T- 
M line, the mirror symmetry forbids the out-of-plane 
spin components and the spin polarization of the sur- 
face states fulfils those expected from the Rashba term 
&R(ky<J x — k x a y ). In contrast, away from the Y-M line, 
the anisotropy of splitting denotes a departure from a 
pure in-plane ordering. The leading correction is due to 
the cubic term (k^. — ky)cr z and results from the in-plane 
deformation of the surface wavefunctions. 

In order to check the spin ordering predicted by the 
measured electronic states, we projected the Kohn-Sham 
wavefunctions on the atomic orbitals of the bismuth 
atoms with defined a z value (z being the surface nor- 
mal). In order to mimic the realistic semi-infinite crystal, 
the summation is focused on the half side of the slab [32 . 
The expectation value of <j y is obtained as <j z but project- 
ing on the spin polarization lying in-plane and perpen- 
dicular to Y-M. Figure |4p and D show the spin polariza- 
tion for the lower and upper surface state, respectively. 
These simulations indicate that the out-of-plane polar- 
ization dominates the in-plane one as soon as <fi > 2°. In 
agreement with our finding, spin-resolved ARPES mea- 
surements reported a large out-of-plane buckling of the 
spins [23] . 

Next, we discuss the entanglement taking place be- 
tween the spin and orbital degrees of freedom near to the 
atomic cores. We obtain the perpendicular component 
of the orbital polarization by projecting of the Kohn- 
Sham wavefunctions on p orbitals with angular momen- 
tum l z = ±1. We verified that contributions from an- 



gular momentum larger than 1 are negligible. As shown 
by Fig. |4p and D, the perpendicular polarization of the 
spin and angular momentum are nearly proportional to 
each other. The scaling factor of -3.5 indicates that spins 
and orbital momenta are counter-aligned. This instance 
resembles the entanglement of J = 1/2 orbitals in the 
strong coupling limit, albeit with electronic states carry- 
ing an orbital polarization quenched with respect to the 
atomic value. Indeed, neither the spin nor the orbital 
momentum is a good quantum number of the system. 
Despite it, the wavefunctions display an average orbital 
momentum L z that directly correlates to S z [24] . 

In conclusion, we investigated the bismuth (111) by 
time resolved photoelectron spectroscopy. The transient 
population of the electronic states above the Fermi level 
allows us to map the SO splitting in reciprocal space. We 
observed a giant anisotropy of the splitting value around 
the high symmetry planes. This finding is not compatible 
with a model based solely on the Rashba coupling and 
chiral spin texture. On the other hand the experimental 
dispersion of the electronic states is in good agreement 
with the results of first principle calculations. It follows 
that a large out-of-plane buckling of the spin and orbital 
polarization is essential to obtain the electronic structure 
of C3 V surfaces with strong SO coupling. Such buckling 
is only a minor perturbation in proximity of the zone 
center but dominates the Rashba term at large electronic 
wave vectors. 
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